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We present a complete set of formulae for calculating the bispectra of CMB temperature and 
polarization anisotropies generated from non-Gaussianity in the vector and tensor mode perturba- 
tions. In the all-sky analysis it is found that the bispectrum formulae for the tensor and vector-mode 
non-Gaussianity formally take complicated forms compared to the scalar mode one because the pho- 
ton transfer functions in the tensor and vector modes depend on the azimuthal angle between the 
direction of the wave number vector of the photon's perturbation and that of the line of sight. We 
demonstrate that flat-sky approximations remove this difficulty because this kind of azimuthal angle 
dependence apparently vanishes in the flat-sky limit. Through the fiat-sky analysis, we also find 
that the vector or tensor bispectrum of B-mode polarization vanishes in the squeezed limit, unless 
the cosmological parity is violated at the nonlinear level. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Recently, the primordial non-Gaussianity of curvature perturbations has been a focus of constant attention all 
over the world. One of the main reasons for attracting so much attention is that meaningful measurement of this 
quantity will become observationally available in the near future, which brings us valuable information about the 
dynamics of inflaton. Bispectrum (three point correlation functions) of the CMB temperature anisotropies has been 
most commonly used to investigate primordial non-Gaussianity [l|, [4] . 

As is well known, if the primordial curvature perturbation deviates from the pure Gaussian statistics, then it 
produces the nonzero bispectrum of the CMB temperature anisotropies. However, there are a lot of sources of 
the bispectrum that not only include the primordial non-Gaussianity of the curvature perturbations but also the 
nonlinearities of the Sachs- Wolfe effect [3- 6] and the radiative transfer l7 Hllll. cosmological recombination G3-QT 
the nonlinear gravitational clustering of dark matter [l6| , the cosmic strings[17l [l8j , the magnetic fields [l9l - l2l| and 
so on. Hence, in order to evaluate the magnitude of the primordial non-Gaussianity of the curvature perturbations 
precisely, it is very important to identify these nonlinear effects. 

In fact, these effects will also induce the tensor and vector-mode perturbations, and the modes may generate more 
characteristic features in the CMB angular spectra than in the scalar one. Some models have been proposed in which 
vector modes are produced at the inflationary phase as well as the scalar and tensor modes, by breaking the conformal 
invariance at that phase [22]]. In such cases, the non-Gaussian vector mode with an interesting amplitude could also 
be generated as in the scalar case. For example, if primordial magnetic fields are considered, the magnetic stresses 
depend quadratically on the primordial Gaussian magnetic field (PMF); hence, their vector or tensor components of 
their bispectra also have finite values. As the vector mode of the CMB transfer function sourced from magnetic fields 
is dominant at small scale, one can expect that the vector-mode bispectrum dominates there in the same manner as 
the power spectra discussed in Refs. [23H25|. Therefore, if one adds the effects of the vector mode in constraining the 
amp litude of PMF by using the CMB bispectra, one will obtain a tighter bound than the current one as O(10)nG 
[20t[2l|. Furthermore, the cosmic strings or the magnetic fields give more characteristic effects also in the polarization 
spectra than in the temperature one |26j-|29[ . Hence, for the identification of the sources of the bispectrum, information 
of temperature and polarization fluctuation generated from tensor and vector-mode ones should be used, not only 
from the scalar mode perturbations. However, there are not enough studies about their effects yet. 

In this paper, we newly present the bispectrum formulae of the CMB temperature and polarization anisotropies 
sourced from non-Gaussianity in the tensor and vector-mode perturbations. First, we formulate all-sky bispectra 
generated from scalar, vector, and tensor modes and find that the bispectrum formulae for vector and tensor modes 
in all-sky analysis formally take complicated forms compared to the scalar mode case due to the dependence of the 
photon transfer functions on the azimuthal angle between the wave vector of photon perturbation k to the unit 
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vector specifying the line of sight direction n. Next, by using the fiat-sky approximation, we simplify the equations of 
bispectra of the CMB anisotropics to solve the above difficulty because no azimuthal dependence arises in this limit. 
In addition, in our flat-sky formulae, we find that if the bispectra of _B-mode polarization is a nonzero value, it infers 
the parity violation in the nonlinear sector. 

This paper is organized as follows: In Sec. [XXI we define the primordial non-Gaussianity from tensor and vector 
perturbations. In Sec IIII| we discuss the formulae of the CMB bispectrum generated from tensor and vector per- 
turbations in the flat-sky analysis. In Sec IIV1 we explain the formulae of their bispectra in flat-sky approximation. 
Finally, in Sec|Vl we provide the summary of this paper. In the Appendices, we derive the formulae of CMB I-point 
function used in the discussion of Sees. Mil and HVl 

Throughout this paper, we assume that the Universe is spatially flat and use the definition of Fourier transformation: 

/(7 3 k - 
(2^ /(k K k ' x , (1) 

/(©) - / ( 2 ) 

where and x are, respectively, 2D and 3D vectors in the configuration space and £ and k are, respectively, their 
Fourier conjugate variables. 

II. PRIMORDIAL POWER SPECTRA AND BISPECTRA OF THE SCALAR, TENSOR AND VECTOR 

PERTURBATIONS 

In this section, we parametrize the primordial non-Gaussianity in the tensor and vector perturbations. As mentioned 
in the introduction, in order to discuss the primordial non-Gaussianity, the bispectrum of the fluctuations is commonly 
used. In this paper, we consider a general expression of the bispectrum of the tensor/vector perturbations which is 
given by 

(r i (k 1 )e(k 2 )e(k 3 )) = (27T) 3 F s ^(k 1 ,k 2: k 3 )S^(k 1 +k 2 +k 3 ) , (3) 

where Sj expresses two helicity states: ±1 for a vector mode, ±2 for a tensor mode. Here, for simplifying numerical 
calculation, we neglect the angular dependence of three wave number vectors in the bispectrum in the bispectrum F. 
This expression includes the so-called "squeezed" or "equilateral" type of the non-Gaussianity [30l l3l| . 

For example, the squeezed-type of the bispectrum is given as follow: As in Refs. [H 0, |32|, [33[ , the primordial power 
spectrum and bispectrum of the scalar curvature perturbations are introduced as 

($ L (ki)$ L (k a )> = (2^) 3 F*(fc 1 ) ( 5( 3 )(k 1 +k 2 ) , (4) 
($ L (k 1 )$ L (k 2 )$ NL (k 3 )) = (2^) 3 F 4) (fc 1 )F$(fc 2 )2/ NL( 5^(ki+k 2 +k3), (5) 

where /nl is the nonlinear parameter of the scalar perturbation and $(k) denotes the Fourier component of the 
primordial curvature perturbation which is decomposed into Gaussian and non-Gaussian part as 

*(x) = $ L (x) + $ NL (x) , (6) 

$ NL (X) EE / NL [$ L (X) 2 -($ L (X) 2 )] , (7) 

in real space. 

This paramctrization can be readily extended to the tensor and vector cases. In contrast to the scalar perturba- 
tion, Fourier modes of tensor and vector perturbations have two independent polarizations. For the convenience of 
calculating the CMB power spectrum as discussed in Refs. [33 - l36| . we use two helicity states (±1 for a vector mode, 
±2 for a tensor mode) to decompose the initial stochastic fields, £ s , where s represents a helicity state. We apply this 
description to the definition of the initial non-Gaussianity of the tensor and vector perturbations as 1 

f(x) = ££(x)+£nl(x) , (8) 
a» = hzUs, -(&'(*)&(*))] • (9) 



1 One can easily include the scalar mode into our notation by considering s = initial stochastic field. For such a case, £° = <!> and 
2-fs 00 = /nl where the index S is used for the scalar mode. 
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where the index Z — T is used for the tensor mode (s,, Sj, Sk = ±2) and = V for the vector mode (s,*, Sj, Sk = ±1). 
Here, we have introduced new nonlinear parameters for the tensor and vector perturbations denoted by fz k s s - • These 
three indices, Sj, Sj, Sk, allow the correlation between each field of the different helicity states in the nonlinear level. 
Because of the symmetry, we have f s z k SiS . — fz k Sj si- By using these expressions, the primordial power spectra and 
bispectra of the tensor and vector perturbations are expressed as 2 

(e(ki)^(k2)) = (2^) 3 ^M<W (3) (ki+k 2 ), 

(e(k 1 )e(k 2 )eN 3 L (k3)) = Pzi 2 k2) /g sis / 3) (ki + k 2 + k 3 ) 

Then, the squeezed type of the non-Gaussianity can be expressed as 

F^( kl ,k 2 ,k 3 ) = ( Pz f l)Pz f 2) /g sl32 +2perms.) . (12) 

In the following discussion, we use the general expression @ as the bispectra of the tensor and vector perturbations 
without specifying the type of the non-Gaussianity. 



(10) 
(11) 



III. CMB BISPECTRUM IN THE ALL-SKY ANALYSIS 



In this section, we derive the formulae of the CMB bispectra sourced from tensor and vector perturbations on the 
full sky. The primordial perturbations in the scalar, vector, and tensor sectors introduced in the previous section are 
transferred through the primordial plasma to the CMB epoch and observed in the CMB temperature and polarization 
fluctuations. In the all-sky analysis, the CMB spin-0 temperature field / and spin-2 polarization fields Q, U, are 
expanded by spin- weighted spherical harmonics [2a . |34| . Following the usual manner, we convert the Q ± iU fields 
into spin-0 E and B fields by using the "spin raising operator" and " spin lowering operator" as Eqs. (|B2[) and (|B3[) . 
Their radiative transfer functions are shown in Appendix [B] 



A. Scalar mode case 



First, we give a brief review of the CMB bispectrum sourced from scalar perturbation. For the scalar case, the 
CMB bispectrum can be written as 

\ a x,i imi a x.e 2m2 a x,i 3 m 3 / - y^Ws °X,l x lzl 3 ' \ La ) 

where the index (5) means that a source of the CMB fluctuation is the scalar perturbation and the index X denotes 
the temperature (/), the -E-mode polarization (E) and the _B-mode polarization (B) 3 . Here, 5™£^ 2 " 13 is the Gaunt 
integral given by 

GiXtr* = J dn y Y eimi (n y )Y e2m2 (n v )Y e3m3 (n v ) , (14) 

and bx £ ± £ 2 £ 3 is the scalar reduced bispectrum formulated as 

F 000 (k u k 2 ,k 3 ) . (15) 

Here T^flih) is the time-integrated transfer function of the scalar perturbation as shown in Eqs. (|B16|) and (|B17[) . 
and ji (x) is the spherical Bessel function. 



b iS) 



V dy 



n 



k'fdkij ti (hy)Tx S J. (K) 



2 For the scalar mode, Ps(k)/2 = P<s>(k). 

3 Of course, scalar perturbation contributes only to the B-mode polarization and not to the B-mode one. In this paper, we use this index 
also for the tensor and vector modes which contribute not only to the E mode but also the B mode. 
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B. Tensor and vector- mode case 



Let us follow the above formulation for the tensor and vector cases. Tensor and vector 1-point functions are 
explicitly given as Eqs. (|B15I) . (|B18I) - (|B23I) in Appendix B. From those equations, one may wonder why tensor and 
vector 1-point functions depend on the spin-weighted spherical harmonics s Yi m (f2fc) although I,E and B modes are 
spin-0 fields. This dependence arises as a consequence of calculating the transfer function in the arbitrary direction of 
the wave vector k. As discussed in detail in Appendix [B] the transfer function for the arbitrary k is written with the 
Wigner D matrix p33l l35l. under the rotational transformation of k from a particular direction (e.g., z direction) 
to an arbitrary direction. This D matrix can be transcribed into s Yl m ((lk) as Eq. (|B13I) . As we will show in the 
following discussion, because of this spin-weighted spherical harmonics s Yi m (Q,k), the CMB bispectra sourced from 
the tensor and vector modes on each angular momentum, £, depends on the sum of the reduced bispectrum over all 
angular momenta 4 . 

For example, let us consider the CMB temperature fluctuation sourced from the tensor perturbation which has the 
spin- 2 spherical harmonics as 

42* => -2^ m (fi fe )£+ 2 (k), +2 Y hn (n k )C 2 (k) . (17) 

Here we consider that the bispectrum of tensor-temperature fluctuations can be sourced from the non-Gaussianity of 
the primordial tensor perturbations which is characterized by the primordial bispectrum given by Eq. ([3]) and hence 
we can easily find 



^LALyS 3m3 ) ^ <5 (3) (ki + k 2 + ks) . (is) 

By making use of the expansion of 3D Dirac delta function given by 



P)(k 1+ k 2 +k3) = / ^e^^)' 



(2tt) 3 

p oo p 3 

= / y 2 dy J dn y J{2Y J i l ' i k-S k iy) Y e^v)Yl'^' i {^k i ) , (19) 
and Eq. (|17j) . a part of the bispectrum of tensor-temperature fluctuations can be expressed as 

i(T) JT) E ^^n/^^miCn*.)-^^*,)- (20) 

ll iLi'„ i=l 



12 3 



Although the bispectrum of the scalar-temperature fluctuation (and .E-mode polarization induced from scalar-type 
perturbation) is derived in the same manner, the orthogonality of the spin-0 spherical harmonic functions gives us 
quite simple expression of Eq. (|13p whose form is the Gaunt integral multiplied by the scalar reduced bispectrum. 
However, as seen in the above expression in the tensor case (also the vector case) the CMB bispectra on each £ depends 
on the sum of the reduced bispectrum over all angular momenta £' as Eq. ([21]) in contrast to the scalar case such 
as Eq. (fT3"|) . due to the nonorthogonality of the 8k dependence between Y^ m i and a Yim (s = ±1 or ±2). One may 
think that this complexity would be evaded once the plane wave could be expanded using spin-weighted spherical 
harmonics, rather than Eq. (|19p . In this paper, however, instead of pursuing this possibility we will use the flat-sky 
approximation to evade this difficulty as we shall show below. 

Thus, the bispectrum formulae of the CMB fluctuations sourced from tensor and vector modes are, respectively, 
given by 



(a {Z) a (Z) a (Z} ) - V /i+4+4c mim2 ™ 3 / v 2 dv 



3 q p oo 

i=l n 



x sgn( Sl ) Sl+a;i sgn( S2 ) S2+ ^sgn( S3 ) S3+ ^F SlS2S3 (fc 1 ,fc 2 ,fc3)^( S ) , (21) 



This complexity does not occur for the CMB 2-point power spectra sourced even from the tensor and vector modes and as is well known 
all CMB power spectra can be described as 

( a jfv'm' a j5m) = c< x>x,e S t't S ™'™ ■ (16) 
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with 

3 . 

ySLi*) = n / ■ ( 22 ) 

i=l 

Here Ty^ is the time- integrated transfer function generated from vector (Z = V) or tensor (Z = T) perturbation 
as described in Eqs. (|B18[) - (|B23|) and x is the index: x = for X — I, E and x = 1 for X = B. 

As we have mentioned before, in the tensor and vector cases, due to the nonorthogonality of the 9k dependence 
between s Ye rn and Y^' m ', the CMB bispectra on each £ depend on the sum of the reduced bispectrum over all angular 
momenta £' as Eq. (j2"Tj) in contrast to the scalar case such as Eq. (IT31 . For this complexity, the numerical calculations 
of the tensor and vector bispectra take much longer time than that of the scalar one. However, this problem can be 
evaded by using the flat-sky approximation as shown in the next section. 



IV. CMB BISPECTRA IN THE FLAT-SKY ANALYSIS 



Here, we explain the formulation of the CMB bispectrum sourced from tensor and vector perturbations by using 
the flat-sky approximation as mentioned in Refs. [f| EH, [34| [Hj]. The flat-sky approximation uses the (2D) plane 
wave expansion of the CMB fluctuation instead of the spherical harmonics one, and it is valid if we restrict observed 
direction n only close to the z axis. As confirmed in Ref. [34f . the flat-sky power spectra of E- and B-mode polarizations 
sourced from the primordial tensor perturbations are in good agreement with the all-sky ones for I > 40. In Ref. [fjj], 
the validity of the flat-sky analysis is also shown in the calculation of the temperature bispectra generated from 
the Sachs- Wolfe term by evaluating the convergence of the modified Bessel function. In addition, in Ref. [l6j], the 
consistency between the flat-sky result and all-sky one in the calculation of the scalar-temperature power spectrum 
and bispectrum are discussed. 

Based on these studies, we have also compared the all-sky power spectra with the flat-sky ones for the /, E, B modes 
from the tensor and vector perturbations and found their consistencies at i > 40. We have also compared all-sky and 
flat-sky temperature bispectra induced from scalar-type perturbations, and confirmed that the flat-sky approximation 
is also applicable in the calculation of the bispectrum for the angular scales where the flat-sky power spectrum is a 
good approximation of the all-sky power spectrum. From these considerations, even if we can not compare the all-sky 
bispectra with the flat-sky ones in the tensor and vector modes due to the difficulties discussed in the previous section, 
we can regard the flat-sky bispectra from the tensor and vector perturbations as good approximations for I > 40. 



A. Scalar bispectra in the flat-sky analysis 

As described in Refs. [H, HI, in the flat-sky approximation the scalar bispectrum Eq. (fT3"|) is modified as 

(a^(e 1 )a^\e 2 )af(e 3 )) = (27r) 2 S (2 \£ 1 + £ 2 + £ 3 )b^ } (h,e 2 ,£ 3 ) , (23) 

Since Q^^ 3 « (2n) 2 S^ {i x + £ 2 +£3), Eq. ^ indicates w 6^(4, 4)- A detailed derivation of 1-point 

functions a^\£) is presented in Appendix C. The scalar reduced bispectra are formulated, by using Eqs. (|C6[) and 
(IC7l). as 



/oo 3 pTq r 00 jl 

y 2 dy J] / d n ^gf{li,h,n,y) 



where Di — tq — Ti and the scalar g functions are described as 



F uuu (fc 1; fc 2 ,fc 3 



gf\l,k,r,y) = S\°>(k,r) 



OS), 



^k 2 - {l/DY D 2 



\A 2 " (e/D) 2 (y - D) 



g { E S) (£,k,T,y) = S£>(k t T) 



- cOS), 



y/k 2 - (I /D) 2 \kDj D 2 



■ cos 



\A 2 - (£/D) 2 (y - D) 



(24) 

(25) 
(26) 



( S) (8) 

Here S T (k,r) and S P (k,r) are the scalar-type source functions of the temperature and polarization fluctuations 
as mentioned in Appendix. [A] 
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B. Tensor and vector bispectra in the flat-sky analysis 



Let us consider the tensor-temperature bispcctrum in the flat-sky analysis. By using Eq. (|C4[) . a component of the 
flat-sky bispectrum of tensor-temperature mode is written as 



(a'-p(£ 1 )aY>(£ 2 )a\ 11 (e 3 )) = 



AT) 



(T), 



n 

.1=1 



dk 



iz a(T). 



OG 



2tt 



1 



-ik ir Di 



(k lz Di) 2 + i\ D\ 



x E pslS2S3 ( v + (v^i) 2 > v fc ^ + (Vd 2 ) 2 , v*^ + (V W) 

Sl,S 2 ,S 3 =±2 



X 



(2vr )3 (5 (2) ( il + fa + fl )(5 (fc l2 + fc 22 + fc 32 ) 

XV 1 1V 2 



By using the expansion of the ID Dirac delta function and the approximation of 2D Dirac delta function as 

S(k lz + k 2z + k 3z ) = ^ ^ e <(*x.+*».+fcs.)» , 

^3) 



* (2) (# + # + #) = D 2 ^ 2 >(£ 1+ £ 2+ £ 3 + ^— _-^ + ---- jD3 



V -Di D 2 D 3 ' ' v Da 

« D?*< 2 >(*i+* 2 + * 3 ) , 



D, 



(27) 



(28) 



(29) 



we can derive the simple form of the tensor-temperature bispectrum in the fiat limit as 

3 



(af ) (£i)af ) (4)af ) (4)) « (2^ 2 )(£ 1 +£ 2 +£ 3 ) / y 2 rfy 

x X F ilMa »(A 1 ,fe2,fc 3 ) . 



TT / / ^- i fl , z T) (4,fci,Ti,2/) 



(30) 



«1,S2,S3 = ±2 



The approximation of Eq. (|29l) is valid because the bispectra are suppressed when the triangle in the ^-space does 
not close as discussed in Ref. 6|. In Eq. (|30|) . we use the approximation (Di/y) 2 1, which is valid because the 
integrand has large value for D\ ~ y ~ ro- Similar to the discussion in the previous section, as the other tensor 
bispectra and the vector bispectra can be derived in the same manner, these bispectra can be written by the same 
form as the scalar bispectra of Eq. (|2"31 which can be written as a 2D Dirac delta function multiplied by the reduced 
bispectra; 



(a%\li)a%\l*)aVp(! 3 )) = (2w) 2 6^(£ 1 +£ 2 +£ 3 )b^(£ 1 ,£ 2) £ 3 ) , 
where the tensor or vector reduced bispectrum is expressed as 



(31) 



y 2 dy 



ft / dr i [ ^r L 9x ) (^ k ii T i>y) 



J2 sgn( Sl fisgn( S2 f 2 sgn( S3 r 3 F s ^(k u k 2 ,k 3 ) 



(32) 



Tensor g functions are written as 



9 P(l,k,r,y) = Sy>(k,r) 



g l P(t,k,T,y) = S^(fc,r) 



v/fc 2 - (£/D) 2 \kDJ D 
kD 



■ cos 



V^ 2 - WW(y - d) 



Vfc 2 - WD) 2 



2 - 



D 2 



■ cos 



^ k 2 _ (t/D fiy - D) 



g<P(£,k,T,y) = -S { P(k,T)— sin [^fc 2 - (£/D) 2 (y - D) 



(33) 
(34) 
(35) 
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Vector g functions are also described as 



sP(k,r) 



cos 



Vfc2 _ (£/D) 2 (y - D) 



i SP(k,T) 



-i SP{k,T) 



kD / D 2 



^k 2 - WW(V - D) 



^(kD) 2 -e 2 d 2 



cos 



\A 2 - WW(v - d) 



(36) 
(37) 
(38) 



Here s\ Z) {k,r) and S K p'{k,T) are the Z-type source functions of the temperature and polarization fluctuations as 
mentioned in Appendix. [X] 

By comparing Eq. (|31l) to Eq. (|23|) , we find that the tensor and vector bispectra are formulated in the same form as 
the scalar one. It is because the helicity dependence, which brings nontrivial couplings between angular momenta in 
the reduced bispectra, vanishes in the CMB 1-point functions induced from the tensor and vector perturbations due 
to the absence of the contribution of azimuthal angle from k to n in the transfer functions as discussed in Appendix [C] 
Hence, unlike the all-sky analysis, the sum of the reduced bispectrum is not needed in calculating the tensor and vector 
bispectra in the flat-sky limit and one can calculate the tensor and vector bispectra with the same computational 
cost taken in the scalar case. This corresponds to the restoration of the orthogonality of Ok between Yp m ' and s Y^ m 
(s = ±1 or ±2) for £ 3> 1, namely, y u i m {s) — > dt^Sz^Sg^ (there is no dependence on mi,77i2 and 7713) due to 
sYim — > Yim for I 3> 1. In other words, it means that because the degeneracy factor of m equal to 2£ + 1 becomes so 
large in the large I limit, the spin eigenstate of s — ±1 or ±2 and that of s — are almost indistinguishable. 

In addition, interestingly, from Eq. (|32[) and (|12j) . we find that £?-mode bispectra (x = 1) from tensor and vector 
perturbations vanish if each fz 3 Sl s 2 ls identical value. This situation corresponds to the parity conservation at the 
nonlinear level. Therefore, if one detects the finite value of the B-modc bispectrum in the squeezed limit, it may offer 
further evidence of the cosmological parity violation. 



V. SUMMARY AND DISCUSSION 



In this paper, we derive the complete set of CMB temperature and polarization bispectra generated from non- 
Gaussianity in the tensor and vector-mode perturbations both in the all and flat-sky analyses. For the primordial 
non-Gaussianity in the tensor and vector sectors, we consider the more general type such as Eq. ([3]), which contains 
the squeezed type given by Eqs. ([8]) and © and the equilateral type. 

Note that the formulation presented can be easily extended in a straightforward manner to the other cases, such 
as a case in which the nonlinear tensor perturbation is excited by the linear-order scalar-tensor couplings. As an 
example of this, we can consider the scalar-graviton interaction during inflation shown by Ref. (39[. Through such 
interaction, the non-Gaussianity of the primordial fluctuations can be generated as a scalar-scalar-tensor type, namely 
) . Although, in the standard slow-roll inflation, such type of non-Gaussianity is expected to be suppressed 
by the slow-roll parameter, it seems interesting that one investigates such type of non-Gaussianity through the future 
CMB observations in the sense of the confirmation of the standard inflation scenario, by using our formulation. 
Furthermore, the 3-point cross correlations between CMB intensity and polarizations, such as (a/a/a^;), and higher- 
order correlations than the 3-point one can be easily formulated in the same manner [33l l40j . 

In the formulation of all-sky bispectra, we find that those formulae take complicated forms compared to the scalar 
one due to the helicity dependence which is represented by the azimuthal angle dependence between the wave vector 
of photon and the unit vector specifying the line of sight direction in the photon propagation. However, in the 
formulation of flat-sky bispectra, we find that the above difficulty is solved for the absence of the above azimuthal 
dependence. In addition, we also show that if the bispectra of B-mode polarization are a nonzero value, it may become 
evidence of the cosmological parity violation in the nonlinear sector. 
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Appendix A: Radiation transfer functions 

Here we show the radiation transfer functions of temperature mode A/ and two polarization modes Aq,Ajj. 
Transfer functions induced by the scalar and tensor modes in a particular basis in which the wave vector of photons 
k is parallel to z axis z are formulated in Refs. [34l l35l . l4lj by the line of sight integral method. For the vector case, 
the method of calculation can be obtained in Ref. [36j Based on the Stokes parameters as defined in Ref. , these 



(Al) 



are expressed as follows: 

AfVo^ftfc.rO - $(k) r dTS^ s \k,T)e-^ 

JO 

A^(r ,k,n fcin ) = (l-/i 2 fe! J$(k) f° drS { p S) (k,r)e-^ x , (A2) 

Jo 

Af)(ro,k,fi fe ,„) = (l-/4 n )(e 2 *^ +2 (k)+ e ^r 2 (k)) / drSf\k,T)e-^ , (A3) 

Jo 

(A<f ± iAP)(t Qi k, n Kn ) = [(1 T ^n) 2 e 2 *-£ +2 (k) + (1 ± Mfc ,„) 2 e- 2 ^-r 2 (k)] 

/■To 

x / dTS { p\k, T )e-^ k ^ x , (A4) 
Jo 

, r T o 

Aj v) (7D,k J fi M ) - -* A /l-Mj I „(f +1 (k)e < **--+r 1 (k)e- < **-) jf drSf ^r)^."* , (A5) 

(Aj ) ± iA<, v) )fo, k, n fc ,„) = ^Hi[T(l T /^.Of"^**'" ± (1 ± mmJ^M^'"] 

x [ ° dTS { /\k,T)e- l >* k ^ x , (A6) 
Jo 

where fifc in (= (Ok,ni 4>k,n)) denotes the orientation of the line of sight direction n in a particular basis in which \c\\z, 

( z) ( z) 

Hk.n = cos 9k_ n . and Sj (fc, r) and Sp (k,r) are the Z-type source functions of the temperature and polarization 



fluctuations [3J, |43 |43J 

In order to estimate the 1-point function ai m , one must construct the transfer functions for the arbitrary k. In 
other words, we want to obtain the transfer functions expressed by the arbitrary k (whose direction is denoted by 
flk) and n (denoted by fl n ) instead of 0/- )n . To achieve this we introduce the rotational matrix 

cos 6k cos (frk — sin <pk sin 6k cos <j)k 
S(Slk)= \ cos ^ sin fe cos0 fc sin 6 k sin cj) k ] , (A7) 
sin 6k cos6>fc 

which expresses the basis rotation that transforms z || k to the arbitrary i. Then the relation between f2fc,n n and 
flk n is written: 



/ sin0 ra cos</)„ \ 




sin0„sin</>„ 


= s(n k ) I 


\ COS<?!)„ / 





(A8) 



COS0t „ 



In the temperature modes, only by changing flk,n to ttk and Q n with the relation (|A8|) . the transfer functions for the 
arbitrary k can be obtained. In the E, B modes, in addition to this treatment, one must consider the mixing between 
Aq and A(j under the transformation S'(fifc) as described in Ref. [Hj]. 
This effect is expressed as 

(A' Q ±iA' v )(r 0> -k,Ct n ) = e TW*(A ±iA ff )(7D,k,n M ) . (A9) 

with the mixing angle tp. The angle "0 represents the rotation angle between 6k, n and 8 n , where 6k, n and 8 n are the 
unit vectors orthogonal to n in a particular basis in which k || z and a general basis, respectively. 

In the flat-sky analysis, i.e., 6 n — > 0, by using Eqs. (jAll) - (|A6[) and (|A8[) and by using the limit of ip as ip — > 
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i>k + 7r, the transfer functions for the arbitrary k are derived as 

Af)(ro,k,O n ) -> *(k) r°dr^ S) (fc,r)e-*-" D , (A10) 

Jo 

(A<f±»A<f ) )(7t )> k I fi B ) -> e* 2i ^"-^sin 2 fe $(k) H drS^ (k,r)e-^ D , (All) 

Jo 

/•To 

Af>fo,k,fi n ) (l- M 2 )(e +2 + r 2 )(k) / dTSi T) (*,r)e- <k - ftD , (A12) 

Jo 

(A^±iA^)(r ,k,ft n ) -> e =F«(*»-*») [(l +M 2 )(e +2 + r 2 )(k)T2M fe (e +2 -r 2 )(k)] 

x / °dT^ T) (fc,T)e- ik -" D , (A13) 
Jo 

/■To 

A { p(T ,k,n n ) i sin0 fe (C +1 +r 1 )(k) / dr5f ) (A;,T)e- k - fiD , (A14) 

Jo 

(a§° ± i A^ ) )(r , k, n„) -> e * 2 ^-^) sin0fe [_ cos e k (t+i + r x )(k) ± - r x )(k)] 

dT^ y) (fc,T)e^ lk -" D . (A15) 



It is important to note that the (f>k dependence which are inherent in the vector and tensor perturbations vanishes in 
the fiat-sky approximation, besides a trivial 4>k dependence due to a spin-2 nature of the Stokes Q and U parameters. 
One may explicitly see that (f>k,n dependence vanishes in the transfer functions when taking 9 n — > because the S 
matrix rotates the basis with the new z axis always being on the x — z plane in a particular basis in which k || z. This 
approximation means that for 5 B <1, it is valid to calculate the CMB fluctuation on the basis of vector and tensor 
perturbations fixed as 9 n — 0, namely, 4>k,n = 



Appendix B: 1-point function in the all-sky analysis 



Here we formulate the all-mode 1-point functions ae m in the all-sky analysis based on the derivation in Ref. [3f 
One-point functions of the /, E, B modes are generated from A/, Aq, Ajj as 



ai,£m — / dVt n 



rf 3 k 
(2tt) 3 



A 7 (r ,k,n n )F; m (fi„) , 



(Bl) 



e?fL, 



(£-2)1 



(£ + 2) 



d 3 k 

1/2 



[(A Q + iAu)(t , k, n n ) 2 Y e * m (n n ) + (A Q - iAu)(T 0) k, n n )_ 2 y/ m (n„)] 



G?fL, 



d 3 k 

(2tt) 3 



8 [A Q + iA u )+d 2 {A Q -iA u ) 



(T ,k,n n )Y; m (n n ) , (B2) 



aB,£m 



~ I dilr, 



(1-2)! 
(£ + 2)1 



d 3 k 

(2tt) 3 

1/2 



[(Aq + iA a )(r 0) k, Q n ) 2 Y; m (Q n ) - (Aq - %A v )(t^ k, fi„)_ 2 y/ m (fi„)] 
d 3 k 



dQ n 



(2tt) 3 



—2 



(Aq +iAu) -d 2 (A Q -iAi 



(r ,k,Q n )Y/ m (Q n ) 



(B3) 



Here we expand with the spin raising (lowering) operators 6 (d ) as introduced in Refs. [U, [13,113 an d ^m(^n) for 
being easily understanding that E,B modes are spin-0 fields, d and 6 act the spin-s function s f{6 n ,<f) n ) as 



d s f(9 n ,4> n ) = -sin s 6>„ [de n +i esc 9 n d<p n ] sin 8 9 ns f(9 n , <j> n ) , 
8 s f(9 n , 4> n ) = -svaT s 9 n [dg n - icsc9 n d rj>n ]sm s 9 ns f{9 n ,(j) n ) . 



(B4) 
(B5) 



Prom here, we derive the 1-point function of tensor-temperature mode as an example. As mentioned in Sec. Mil 
This is calculated by using Wigner D- matrix Di^ , , which is the unitary irreducible matrix of rank 2£+l that forms 
a representation of the rotational group. The property of this matrix and the relation with spin- weighted spherical 
harmonics are explained in Refs. (33l. [35l \3l\. By using Eq. (|B1|) . the relation between the Yt m and D matrix, and 
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the relation corresponding to Eq. (|A8I) as 

*7m(n») = J2 D( rnL(s(n k ))Y; m ,(n k , r 

m' 

the 1-point function of tensor-temperature mode is written as 

d 3 k 



>I,tm- J (27r)3 

Next, with the mathematical relations as 

1 d/i fc , n (l - ^.n)^ 2 (^,n)e-^'" X 
-1 

the integration for £7fc in can be performed to obtain 



2/ + l(l-m')! 
4tt (£ + m')! 
(£-2)! 



1/2 



*T (/**,») e ~ 



(£ + 2) 
-2(-i) 



-Pjf (/*fe,n) j 

,(£ + 2)!^) 
(i-2)! ' 



(B6) 
(B7) 

(B8) 

(B9) 
(BIO) 
(Bll) 



A T ) 



~{l + 2)\' 


1/2 


"2^ + 1" 


1/2 /■ d 3 k 


L(^-2)!. 




47T 


J (2tt)3 



^2 (5(n fc )) e +2 (k) + i>w ( S (n fc )) r 2 (k) 



/o ^ 



Because £) matrix is written by the spin-weighted spherical harmonics as 



47T 

2£+ 1 



1/2 



(B12) 



(B13) 



we obtain the final form, namely 



AT) 

l IJm 



'{£ + 2)1' 


"7 


W-2)\_ 





(27T) 



[_ 2 y; m (^.)e +2 (k)+2^„(O fc )r 2 (k)] H drSf \k,r) 3 -^- . (B14) 



For the other modes, we can derive in the same manner with Eqs. fBTj) - (|B3)) , (|A1|) . ([X2]) . (|A"3]l - (|A"6)l , 
As a result, all-sky 1-point functions can be formulated: 



'Xim 



4 7 r(-i) < 



d 3 k ^ /^(fi*)*^)^?^) 



(for Z = S) 



(2*) 3 lE s sgn( S )^_ s Y/ m (n fc )e(k)7l7(fc) (forZ = T,V) ' 



(B15) 
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where x = for X — /, E, x — 1 for X = B, time-integrated transfer functions Tp){k) are expressed as 



7#>(fc) = / drSP(k,rMx), 



r<V(k) 



Tfhk) 



■2V. 



(1 + 2)1 



1/2 pro 



-2Y. 



{i-2)\ 

To 



" drS { p^ s \x)j t {x) 



V2 r d rSP(k,r)^ 



o 



drS^(k,r)£^(x)^ 



dTS { p(k,T)B (T) (x) — 

ar 



[(^+1)!" 


1/2 


[(£-!)!_ 


J 


(£ + 1)! (£ 


-2)1 


(£-l)\(£ 


+ 2)! 


(£ + l)\ {I 


-2)! 




+ 2)! 



r ° drS\ V \k,r)^l 



1/2 



1/2 



' d rS { P{k,r)£^(x) ] ^ 
x 

'drSP(k,r)B^(x)^ 



(B16) 

(B17) 

(B18) 
(B19) 
(B20) 

(B21) 

(B22) 

(B23) 



and the operators £, B are defined as 

£^(x) = (l + d 2 x ) 2 x 2 , 

£^{x) = -12 + x 2 (l-8 2 x )-8xd x , 

B (T) (x) = 8x + 2x 2 d x , (B24) 
£^ (x) =Ax + (12 + x 2 )d x + 8xd 2 x + x 2 d 3 x , 
B {v \x) = x 2 + Axd x + x 2 d 2 . 

Note that in the all-sky analysis, due to the dependence of transfer functions on 4>k,m 1-point functions depend on 
the helicity state through the spin spherical harmonics. 



Appendix C: 1-point function in the flat-sky analysis 



In this section, we formulate the all-mode 1-point functions ag m in the flat-sky analysis. In this limit, 1-point 
functions in the all-sky analysis described as Eqs. (|Blj) - (|B3P are modified by using the plane wave as 



Q-llra 



^2 ld2& 



— > 



d 2 & 



d 3 k 
(2nf 
d 3 k 
(2tt) 3 
d 3 k 
(2ir) 3 



(A Q + zAu)e' 2 ^-^ + (Aq - iAu)e 2 ^-^} (r ,k, fi„)e-* = a E (£) , (C2) 



-(A Q +iAu)e 



(A Q -iAu)e 



2i(4> e -4> n ) 



(T Q ,k,n n )e- l£ - @ =a B (£) ,(C3) 



where is the 2D vector projecting n to the flat-sky plane expressed as © = (0cos0„, 0sin</>„). For example, in 
order to obtain the 1-point function of the tensor-temperature mode, we substitute Eq. (|A12[) into Eq. (IC1|) and 
calculate as follows: 



d 3 k 
(2^)3 

d 3 k 
J2nf 



a ^ [i) = J J^^-tiH^+C 2 )^) ['dr J d 2 ®e-^ lD +*y & S P(k,T)e-^ D 
i 2 k (C+ 2 + C 2 ) (k) f dT(2n) 2 S^(^D + t)sf\k, T)e~^ D 



/•To /-OO 77 p2 1 

I dT J x fr (e+2 +r2)(k " = - l,D ^ \k z Df+p Si (fc = ^+W^)^ e_,lifl '( C4 ) 
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where D = To — r is the conformal distance and we have decomposed k into two-dimensional vector parallel to the 
flat sky and that orthogonal to it, k = (k", k z ). In order to obtain the last equation, we use following relations which 
are satisfied under k" = —£/D as 



k = \lkl+ ( - 



£ ,2 



sin V k 



J_ _ i 

kD ~ ^{k z DY+P ' (C5) 



cos6» fe = sgn(fc z )y 1 

4>k = $>l + 7T . 



2 

kD . 



One-point functions of the other modes are calculated in the same manner by using Eqs. (|A10I) . (| Al 1|) . (| A13I) - (IA15[) 
and EH) - 0) as 



PTq poo JI, -| 

a\ S) (£) = j drj ^$(kH = -£/D,k z )Sf\k = y/k* + {l/D)\r) — e^ D , (C6) 

f T ° f 00 rlk P 2 , cl , 1 

l * L ^ $(k = = V^(W^)- e -- , (C7) 



af{t) 



TO 



a%\£) = /"dr/ ^(^+r 2 )(kH=-Wfc z ) 



a ( P(£) 



x ( 2 - (kjy+p ) ^ T)(fc = v^FT(W, , (C8) 

pro roo 

i dr dk z (t 2 - r 2 )(k !l = -l/D, k z ) 

JO J-oo 

x2 sg n(k z ) ] jl^^^S { P(k = v/fc, 2 + (^) 2 ,r)-^ e -^ , (C9) 

/•To r oo 7 7 

4 v) w = * / dr ^(e 1 +r 1 )(k" = -//Afe) 

JO J-oo ^ 

S [ P{k = ^/kJTWW,r)^-e-^ D , (CIO) 



«2°w 



/•TO /-OO 77 

/ dr/ p(e +1 +r 1 )(kN = -^/A*,) 



x S gn(^)- 7 =^_^/l - j^^S^ik = VWTWW,r)±e-^ , (Oil) 



TO 



a<P{£) = -i J "dr j°° ^(^ -r X )(kH = ~£/D,k z ) 

x - £ =S ( p \k = ^kl + {e/D)\r)^e- tk * D . (C12) 

Note that the helicity dependence vanishes in flat-sky 1-point functions unlike in the all-sky ones as shown in Appendix 
151 It is due to the absence of 4>k,n dependence in the flat-sky transfer functions as explained in Appendix lAl 
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